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Screening of topological charges (singularities) is discussed for paraxial optical fields with short 
and with long range correlations. For short range screening the charge variance (Q^) in a circular 
region A with radius 7? grows linearly with R, instead of with R^ as expected in the absence of 
screening; for long range screening (Q^) grows faster than R: for a field whose autocorrelation 
function is the zero order Bessel function Jo, {Q^) ^ RlnR. A Jo correlation function is not 
attainable in practice, but we show how to generate an optical field whose correlation function 
closely approximates this form; screening in such a field is well described by our theoretical results 
for Jo. (Q^) can be measured by counting positive and negative singularities inside A, or more 
easily by counting signed zero crossings on the perimeter P of A. For the first method (Q^) is 
calculated by integration over the charge correlation function C (r), for the second by integration 
over the zero crossing correlation function F (r). Using the explicit forms of C (r) and of F (r) we 
show that both methods of calculation yield the same result. We show that for short range screening 
the zero crossings can be counted along a straight line whose length equals P, but that for long 
range screening this simplification no longer holds. We also show that for realizable optical fields, 
for sufficiently small R, (Q^) ~ R? , whereas for sufficiently large 7?, (Q^) ~ R- These universal 
laws are applicable to both short and pseudo-long range correlation functions. 



I. INTRODUCTION 



Random (and other) paraxial optical fields generically contain numerous point topological singularities (defects) in a 
plane (the a;?/-plane) oriented perpendicular to the propagation direction (2;-axis). These singularities, which include 
phase singularities (optical vortices) [TJ [5] polarization singularities (C points) [3] , gradient singularities (maxima, 
minima, and saddle points) [4], and curvature singularities (umbilic points) [5], are the defining features of the optical 
field, and are characterized by signed winding numbers (topological charges) q±. Generically, vortices and gradient 
singularities have charge q± = ±1, whereas the charge of C points and umbilic points is q± = ±1/2. 

Like electrostatic charges, topological charges screen one another: positive charges tend to be surrounded by a net 
excess of negative charge, and vice versa [6l III [3 [9l [lOl [TTJ [12j [131 IHl US [El Ull IE] ■ A formal measure of screening 
is the charge correlation function C(r) [B], which by convention is constructed to measure the net excess of negative 
charge surrounding a positive charge at the origin. When 

j j C{r)dr — —1 (^ 2 ^ ^^'^ umbilic points^ (1) 

screening is said to be complete, whereas for partial screening — 1 < JJ^^C{r)dr < 0. Here we concern ourselves 
with complete screening. 

In generic optical fields the average net charge (Q) = , but there are large fluctuations: these are characterized 
in lowest order by their variance (Q^)- For a purely random collection of N singularities with say charges q± = ±1, 
(Q^) ^ N. Screening significantly reduces these fluctuations, to an extent which depends on the rate of fall-off of C(r) 
for large r. In the generic case where C(r) decays sufficiently rapidly, canonically exponentially, (Q^) ~ N^^"^ [S]; this 
is the hallmark of short range screening,. For a slowly decaying C(r) however, we shall see that the charge variance 
grows faster than N^/'^; we call this regime long range screening. 

(Q^) and C(r) are closely related, and given the latter the former can, in principle, be calculated. Here, we 
explore in detail the behavior of (Q^) and its relationship to C(r) for long range screening. Speciflcally, we examine 
a particular example, first suggested in [10 , and obtain analytical results for the behavior of (Q^) in a large circular 
region of radius R, finding (Q^) ~ RlnR. 

In elliptically polarized paraxial fields C points are vortices (zeros) of the right and left handed circularly polarized 
components of the field [3], so our results are applicable also to these singularities. C{r) is also known for the 
stationary and the umbilic points of random circular Gaussian fields, such as the real and imaginary parts of the 
optical field [15], but the correlation functions of these singularities differ from that of the vortices and our results are 
not necessarily applicable to these singularities. 

Measurement of C(r) (either from experiment or computer simulation) requires locating and characterizing all 
vortices in each of a large number of independent realizations - this is a challenging task. On the other hand, Q, 
and therefore Q^, can, using the index theorem be obtained from measurements made only on the boundary of 



2 



the region of interest. As discussed in [S], a convenient way of accomplishing this is to count signed zero crossings of 
either the real or imaginary parts of the optical field. Correlations between these zero crossings are described by the 
zero crossing correlation function T (r). In [3] it was argued that like C(r), also T (r) can be used to calculate (Q^)', 
this argument was supported by numerical results and computer simulations, but was not demonstrated analytically. 

The plan of this paper is as follows: In Section II we discuss the connection between (Q^) and C(r) in a natural, i.e. 
sharply bounded, region of finite area - first reviewing known results for short range screening, and then presenting 
what are to our knowledge the first explicit results for such boundaries for long range screening ^. These latter are 
the major contribution of this report. In Section III we discuss the connection between (Q^) and the zero crossing 
correlation function r(r), showing, analytically for the first time, the equivalence of the two seemingly different 
approaches to {Q^}'- one based on C(r); the other on T (r). In Section IV we consider a practical realization of an 
optical field with long range correlations that could be used to compare experiments with the results presented here. 
In the discussion in Section V we consider asymptotic forms of {Q^} for very small, and very large, R, for physically 
real optical fields. We summarize our main findings in the concluding Section VI. 



II. CHARGE VARIANCE IN A BOUNDED REGION AND THE CHARGE CORRELATION 

FUNCTION C (r) 

In [5] the following quantitative relationship between the charge variance (Q^) and the charge correlation function 
C(r) was derived for a bounded region of area A Eq. (37)) with ({ANf) = (Q2)]: 



{Q^) = {N) + ^ ir A{v)C{v)dv. (2) 

Here and throughout (...) represents an ensemble average, N is the number of singularities with charges ±1 contained 
in A, rj is the average number density of singularities, and A{r) is the area of overlap between A and its replica 
displaced by r. 

In what follows we assume isotropy, and therefore take A to be a circular area with radius R. Elementary geometry 
yields 

A{r) ^ A{r) = TrR^ + B{r), (3a) 



Bir) . -2i?^(^arcsin(^) + ^^l-(^)^j, 0<r<2R, 

= 0, r> 2R. (3b) 
Noting that (N) = TyTri?^, Eq. ^ becomes 





p2R 






l + 27r / rC{r)dr 


+ 2717] [ 




Jo 


Jo 



r2R 



(4) 



Eq. ^ is our basic starting point. We proceed to evaluate it for large R: first for a characteristic field with short 
range correlations; then for one whose correlations are long range. 



A. Charge correlation function 

Halperin [B] first showed that assuming stationarity and circular Gaussian statistics, the charge correlation function 
for optical vortices can be written as 

^ W'jr) [W{r){W'{r))^ + W"{r){l - W'jr))] 
7rM^"(0)r(l-M^2(^))2 ' 



A long range screening result for a Gauss-smoothed boundary is given in [lOj . 
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where 



W{r) = {EXO)E{r))/(\Emn, 



(6) 

is the normahzed autocorrelation function of the optical field E, and W'{r) — dW{r)/dr, W"{r) = <PW {r) / dr"^ . Liu 
and Mazenko |7j, and Berry and Dennis [10 , subsequently noted that Eq. ([5]) can be cast into the convenient form 



C{r) 



where 



-1 dO? (r) 
2-KrW"{Q) dr ' 

W'{r) 



^1 - W^{r)' 

The number density of vortices 77 also depends on W{r), and is [2] 

~-W"{Q) 



27r 



A useful measure of how rapidly screening sets in is 



Ic{p) = 2ii / rC{r)dr 



W"{Q) 



1 - W^{p) 



(7) 
(8) 

(9) 

(10a) 
(10b) 



For short range screening / w — 1 for p of order nearest neighbor separations, or greater; for long range screening Ic 
approaches —1 only when p — > 00. 



B. Short range screening 

If C{r) decays sufficiently rapidly with r, the large r contribution in the first integral in Eq. Q is negligible, and 
the upper limit can be extended to infinity. Using the screening relationship in Eq. ([T]) we then have 

(Q2) ^ 27r77 / rB{r)C{r)dr. (11) 



In this same limit B (r) in Eq. ( 3b I can be replaced by its leading term in an expansion in powers of r / R, B{r) —2rR, 
and we have 

(q2) = l^A.P, (12) 

where P — 2ttR is the length of the perimeter of A, and the screening length As, which is a measure of the width of 
Icip) inEqs. (10), is [Hj 



As^-S / r^C{r)dr. (13) 

Thus, the hallmark of short range screening is that for R >> A^, (Q^) grows with the perimeter, rather than with 
the area, of A. 

These results have a simple physical explanation [S]. Because charges deep inside A are perfectly screened they 
make no contribution to the net charge Q, and therefore no contribution to its variance {Q^}- Charges located within 
a distance of order A^ from the boundary of A, however, are imperfectly screened. The number of these charges 
is n « rjAgP. For completely unscreened charges (Q^) ~ n, whereas for partially screened charges we can expect 
(Q^) « an [S], where from Eq. ([12]) a= \. 

Eqs. ( [T2| and ([l3| also quantify the notion of short-range screening. The condition that A^. is finite is sufficient 
to insure that {Q'') grows linearly with R. Thus, for short range screening C{r) (averaged over a period in the case 
that it oscillates) must fall faster than 1/r^. 
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1. Gaussian correlations 



The canonical short range correlation function which gives rise to short range screening is the Gaussian, which we 
write as 

W{v) = exp (-K^a^w^) , (14a) 
K = 27r/(AZ), (14b) 

where A is the wavelength, Z is the (asymptotically large) distance between the random source and the screen on 
which the field E is measured, and v measures radial displacements on this screen. 

W{v) is the normalized autocorrelation function of the far field speckle pattern produced by a distribution of 
randomly phased sources with amplitudes 

5(M)=exp(-K(2a)]'), (15) 

where u measures radial displacements in the source plane, and a is a measure of the width of the distribution ^ . The 
displacement r and the displacement v are related by the scale factor k 



giving for the Gaussian 



Then, using Eq. ( 13 1 we obtain 



r = /cv, (16) 
1¥ (r) = exp (-aV^) . (17) 

= C (3/2) / [V2^a^ = 1.0422/a, (18) 



where C is the Riemann zeta function. 

The average spacing d between singularities is 

d= ^/lfq = ^fn/a=l.n/a, (19) 

which exceeds the screening length Ag by nearly a factor of two! This apparent paradox is resolved in Fig. [T| where 
it can be seen that the singularities tend to cluster with nearest neighbor spacings of order A^. 

A measure of how quickly screening sets in, thereby leading to local charge neutrality, is the rate of decay of 
Ic (p) + 1, Eqs. (10): for p = A^, /c + 1 = 0.3, for p = 2As, it equals 0.003, and for p = 3As, /c + 1 has dropped 
below 3 X 10"^. 



2. Smoothed boundaries 



In ([To], Eq. (4.^ 



for the charge variance, here (Q^) 



and (4.49)) a Gauss-smoothed boundary is assumed, and a constant, i.e. R independent, value 
is obtained from 



N 
~2 
1 



1 + 2717] / X exp 







TTX 



gQ{x)di 



4./o l-W^{x) ' 



(20a) 
(20b) 



where C (r) in [10 is our W (r), gq (x) — C [x) /rj with C [x) given in Eq. (|5]), A is the area of the region of interest, 
here A = 7ri?^, and N = ijA is the total number of charges inside this region. Eq. (20b I is obtained from Eq. (20a I 



^ Such a source distribution can be obtained in practice by illuminating a random phase screen (typically a piece of finely ground glass) 
with a Gaussian laser beam. 





(a) 



(b) 



FIG. 1: Singularity structure. Shown is the phase of the optical field coded to 2n black to white for a field with (a) a 
Gaussian and (b) a Jo correlation function. Positive (negative) singularities are shown by filled white (black) circles. Both 
fields have the same average number density of singularities, (a) Gaussian correlation function. Here there are 13 positive 
and 13 negative singularities. The screening length As is shown by the bar. Note the clustering of singularities into small, 
charge-neutral groups with the spacing between nearest neighbor singularities being of order As, and with large empty regions 
between groups. It is this clustering that produces a screening length ~ 1/2 the average spacing between singularities, (b) Jo 
correlation function. Here there are 14 positive and 14 negative singularities. As can be seen, there is a strong tendency for the 
singularities to be nearly equally spaced on a lattice-like structure in which each nearly square unit cell, shown by the dotted 
rectangle, contains two positive and two negative singularities. This structure exhibits substantial local charge neutralization, 
but because of defects in the lattice the overall screening is long range. 



by using Eqs. ([T]), ([t]), (8), and (joj), expanding exp {—ttx^/ (2/1)) keeping the leading term in x, and integrating by 
parts. 

For the Gaussian correlation function in Eq. ( [T7| , Eq. ( |20b| can be evaluated analytically, and we obtain 



0.2056. 



(21) 



This result shows that (Q^)j,o„sj is not only independent of R, but also of a, and is therefore independent of the charge 
density r] and all other system parameters! A question of considerable interest not discussed in [TU] is whether or not 
this value for (Q^)^^^^^ is an intrinsic property of the medium, or is it the result of a particular set of assumptions? 



In [To] it is stated that Eq. (20a I is obtained from a calculation "with the boundary of A Gauss-smoothed to 



eliminate trivial edge effects" , but the details of the calculation are not given. In order to be able to more closely 
examine the physical content of Eq. (20a) we need a derivation of this equation. Here we briefly outline such a 



derivation, starting with the deflnition of a boundary that is "Gauss-smoothed" . 
A{x) in Eq. (3) can be obtained from an area function a{v) using 



(x) = 2 / rdr I dO [a (r) a (r -I- x)] 
Jo Jo 



(22) 



Throughout this report we have used a disk shaped area described by the area function 

a(r) = a(r) = 1, < r < R, 
= 0, r>R. 

This function, shown in Fig. |2]as curve (a), inserted into Eq. (22 1 yields 

/>arccos(2;/(2r} 

^ (a;) = 4 / rdr / d9, 

Jx/2 Jo 



(23) 



(24) 



which when evaluated yields Eq. (3). 
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FIG. 2: Area functi ons , 
area function, Eq. (23 1, 
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p = 0, a so called super Gaussian; curve c, p 
so called stretched exponential. 



Shown is the normalized area function a{r) for different values of p in Eq. ( |29[ |. Curve a is the disk 
formally Eq. (291 with p oo. Curves b - e are Eq. (29 1 with the following values for J : curve b, 

1/2, a 



2, a normal Gaussian; curve d, p = 1, a normal exponential; curve e, p ■ 



We note that in order to maintain the same number of charges N within A as there are for the disk, any other form 
for a (r) must satisfy 



2tt i ra (r) dr = A = ttR^ 
Jo 



'0 

For a Gauss-smoothed boundary we have instead of Eq. ( 24 1 , 

„2 , 



a 



(r) = exp (-r^/i?^) = exp [-nr'^/A) 



which is shown in Fig. 2 as curve (c). Inserting Eq. (26 1 into Eq. (|22[) yields 



2A 



A{x) = -Aexp - 



(25) 



(26) 



(27) 



A peculiarity of Eq. (27) is that -4(0) = A/2, as opposed to -4(0) = A, the expected result, which is obtained 
from Eq. (3). As a result of this peculiarity, in the absence of screening the Gauss-smoothed calculation yields 
{Q'^) = N/2, instead of the expected random walk result (Q^) — N. 

Generalizing Eq. ^ to an arbitrary area function satisfying Eq. (251, we have 



(Q2) = r]A (0) + 2Tn] / xA {x) C {x) dx 



(28) 



Inserting Eq. ( 27 ) , and using C (x) = rjgq (x) we recover Eq. ( 20a I . 

Now, if the physical content of {Q'^}^gnst '^^^^ that it is an intrinsic property of the random medium, rather than 
simply the result of a calculation for a particular form for a (r), then other forms for a (r) should yield the same result. 
In order to test this we redo the above calculation using a boundary smoothed by a general exponential function of 
the form, 



a (r) = exp {-bp {r/Rf) , 



where < p < oo is a real, not necessarily integer, number, and in order to satisfy Eq. (25 1, 



bp = 















p 







(29) 



(30) 
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Examples of a ir) for different values of p are shown in Fig. 2 as curves (b), (d), and (e). 
Inserting Eqs. (|29| and Q into (|22|, we find 



(31) 



so that in the absence of screening (Q^) = 4^(^/P^iV, which approaches N for large and approaches zero as p 
approaches zero. 

For the charge variance itself we have 



(32) 



which reduces to the Gaussian result in Eq. (21 1 for p = 2, but differs from it for other values of p, demonstrating 
that (Q^)^^, is, in fact, not an intrinsic property of the medium, but is dependent on the particular choice of a (r). 
In [3] a short range autocorrelation function of the form 

W{r) = 2Ji(ar)/(ar) 

was used. For this form for arbitrary p we obtain by numerical integration 

(Q')eo„.*-P (0.168077...), 

which is again independent of a. 



(33) 



(34) 



Because both Eqs. (32) and (34l are continuous functions of p, a particular value for (Q^)^ 



, , does not uniquely 

' const ^ 

specify neither the field autocorrelation function nor the boundary smoothing function. For example, {Q"^) ^^^st ~ 
0.2056... for a Gaussian autocorrelation function and a Gaussian boundary smoothing function for which p — 2. But 
the exact same value for is obtained using the Ji autocorrelation function in Eqs. (33 1 and a boundary 

smoothing function with p = 1.1817.... 



C. Long Range Screening 



Berry and Dennis |10j have introduced the long range correlation function 

Wir)^Joiar), 



(35) 



where J„ is a Bessel function of integer order n. This form for W (r) arises from a random source distribution S (u) 
that takes the form of a very thin ring. 



5* (u) = (5 (u — a) . 



(36) 



Practical realization of a Jq correlation function is discussed in Section IV. 

Using Eqs. (10) with p ^ oo the screening relationship in Eq. ([T|) is easily verified, whereas upon inserting Eqs. 
([t]) and (|8| into Eq. ( [T3| it is immediately seen that As does not converge. These are the conditions for long range 
screening. 

Inserting Eqs. ( |3b[ ), ([T]), and ([s]) into Eq. Q, and integrating by parts, we obtain 



27r Jq 1 



{ W'{r)f 
W2(r) 



dr 



(37) 



In the remainder of this report we use this highly convenient form for {Q"^) exclusively. 

For the case of Jq, Eq. (|35|), we break the integral on the R.H.S. of Eq. (37 1 into two parts: 



3l +^2, 



where 



3i = 



32 



a2 


f\ 


27r 


Jo 


a2 


I-2B. 


2tt 


J A 



(38) 

(39a) 
(39b) 
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with 1/a < A < R. 

In evaluating 3i and J2 we will need the asymptotic forms 

Jiiar) o, , 2 cos^ (ar + 7r/4) _ 
. ,2/ ^ ~ -^1 K ~ ^ r » 1/a. 

Expanding \/4i?^ — r"^ and keeping only the leading R term we have for Ji, 



ai? , Jl{x) 
/ "2; „ 

This form diverges as a A — > 00, and so we regularize it by writing 

Jl{x) 



3i 



aR 






1 dx 


T J 











TT 


Jo 



1 - J2(:r) 



aA 



where 



' ^oix)JH^) ^ 0.5630468586... 

i-'^o(a;) 



(40) 
(41) 

(42a) 
(42b) 

(43) 



is evaluated numerically. Using Eq. (40 1 the remaining integral in Eq. (42b I can be evaluated analytically in the 
limit of large a A as follows: 



aA /-aA 

J?(x)dx — lim / dx x~'^J?(x), 

£^0+ ./n 



lim 



lim 

0+ 



lim 

£^0+ 



dxx ^Jf{x)- / dxx "Jfix) 

JaA 



2FI 



2^r(2)r(i±s) 
r(5=^) r(2)r(6) 



^,^;2;l 
2 ' 2 ' ' 



qA t^x 



2£r(2)r(i±^)r(i±^)r(^) 

7 + 3 In 2 - 2 1 



(aA)- 



In(aA), 



(44a) 
(44b) 

(44c) 

(44d) 
(44e) 



where 7 = 0.5772... is Euler's constant. 



Similarly, in calculating Eq. p9b|, we use the asymptotic form in Eq. (|40|) to obtain 

dr 



J2 



27r2 
a 

a 
27r2 
aR 



2« V4i?2 - r2 



A 



V 4i?2 - a;2 - 2i? In 



a: 

8i?2i 



2i? 



A 



4i?2 

-2i?ln 2i? + 2i?ln , 

2i? A 



A 



Combining 3i and ^2, we obtain the central result of this report, 



aR 



= — [K + \n{aR)] 



(45a) 

(45b) 

(45c) 
(45d) 

(46a) 
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where 

if = 7rI? + 7 + 51n2-3 = 2.81181544... (46b) 

We note that for rather modest values aR > 17 the In (ai?) term dominates, so that this term should be easily 
accessible to experiment. Direct numerical integration of Eq. ^ reveals that Eq. (46a) is good to better than 3% 



for aR = 2, to better than 0.3% for aR = 10, and that thereafter the percentage error decreases as '-^ 3/{aR). 

As for the Gaussian, a measure of how quickly (or slowly) screening sets in is the average rate of decay of Ic (p) + 1. 
For large p this average rate is ~ 2/ (irap), showing that although screening is ultimately complete, it approaches 
completeness very slowly. 



1. Smoothed boundanes 



What does boundary smoothing do to (Q^) when the correlations are long ranged? 

For the long range correlation function Jo(ar), [10] states "we can show from (4.48) [here Eq. (20al] that . 
(Q2) ^ y/N", i.e. for a circular region of radius R, {Q"^) ~ R. 

For large i?, the limit also used in [10 , we find here for all generalized exponential smoothing functions in Eq. ([29 



{Q^)^^{p)aR. (47) 
For the Gaussian smoothing function used in [10] we obtain analytically 

S (2) = = 0.09974. (48) 

4v27r 

In general, however, S (p) must be evaluated numerically. Introducing scaled variables x = hy^x/R,r = by^r/R, 
and using the asymptotic form fi^ [x) ~1/ (ttx) , Eq. ( 40 ) , we have 



2TTHy" Jo \ dx , 

A{x) = 2 ¥dr d6lexp(-7~''-fr^), (49b) 
Jo Jo 

f± — + x^/4 ± rx cos 6. (49c) 

Carrying out the required numerical integrations we obtain: S (4) = 0.1615; S (6) ~ 0.2052; S (8) = 0.2338; and 
S(10) = 0.2561. 

Thus, also for long range correlations the results of boundary smoothing change when the arbitrary smoothing 
function is changed. 



III. CHARGE VARIANCE IN A BOUNDED REGION AND THE ZERO CROSSING CORRELATION 

FUNCTION r 

As indicated in the Introduction, Q, and therefore (Q^), can be obtained from measurements made only on the 
boundary of the region of interest A, here a circle of radius R. The method suggested in [5], which is amenable to 
calculation, is the following: Signed zero crossings (ZCs) of either the real TZ, or of the imaginary X, parts of the 
wave functions that cross the boundary are counted, with each positive (negative) ZC contributing +i (—5) to Q; 
the sign a-jz of a zero crossing of TZ is a-ji = sign(— XT^^), that of a ZC of 1 is aj — sign{TZXs), where 71^ = dTZ/ds, 

Is = di/ds [g. 

Writing the lineal number density of positive (negative) zero crossings of 7?, or of X as n+,(n_), their sum and 
difference are 

no = (50a) 
AiV = n+-n^. (50b) 

where, for circular Gaussian statistics [T^] 



no= -^y-W" (0). 

TT 



(51) 
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For a straight line of length L oriented along say the x-axis [S], 



The zero crossing correlation function F is 



F(Aa;)fi(Ax) 



L 



|Aa;|F(Aa;)d(Aa;) 



1 (1- VF2(Aa;)lF"(Ax) + W(Aa;)(VF'(Ax))2 . , 
dfl{Ax) 



^^^^^ TT^no (1- VF2(Aa;))3/2 



TT^no 



arcsin (W (Aa;)) 



d{Ax) ' 



(52) 

(53a) 
(53b) 



where is obtained from Eq. (|8| with r replaced by Ax. F(Aa;) and C (r) are compared in Fig. [sj 
For a circle of radius R, 



1 



no (27ri?) 



F(A6')d(A6') 



(54) 



where A9 measures the angular separation of two points on the rim of the circle. This form for ^(AiV)^y follows 
from Section F of [9 with x replaced by 0. The corresponding form for F (AO) is 

F(A6') 



2 W 



1 . dO(A0) 

arcsm [W) 



d[AB) ' 



(55) 



wher e tt-q is the number of zero crossings per unit of arc that cross the rim of the circle, and il{A9) is obtained from 
Eq. (53b I with Ax replaced by A9. This form for F {A9) follows from Section G of jS] with x replaced by 9. 

Using the fact that in the limit A9 goes to zero the difference between the rim (arc) of the circle RA9 and the 
corresponding chord, here r, vanishes, together with Uq \d9\ = ng jdrj, we have 



A0) 



Rno, 



(56) 



where no, given in Eq. (50a I, is the lineal number density of zero crossings. 



In what follows we refer to Eq. (52) as the linear F formula, and to Eq. (54 1 as the circular F formula. 



The field autocorrelation function W {A9) appearing in ^{A9) depends only on the straight line separation, i.e. 
chord length r, between two points on the circle rim separated by A9, 



W{A9) = W{r), 

r = 2Rsm{A9/2) 



(57a) 
(57b) 



A. Exact Equivalence 

We now proceed to show explicitly that (Q^) in Eq. [2j which is written in terms of the charge correlation function 
C and the area of A, and (Q^) in the circular F formula, Eq. (54 1, which is written in terms of the zero crossing 
correlation function F and the perimeter of A, are identical - as they must be. 

Returning to Eq. ( 37 1 we write 



(0") = ^ / dr V4i?2-r2 — 
27r Jo Wl 



W'{r) 



^1 - VF2(r) dr 



1 

2^ 



TT RUo 



2R 



dr sin ^ W{r) 



dr 



sin-^ W{r), 



VFR2 



W'{r) 



^1 - VF2(r) 



We now compare this with the circular F formula, Eq. (54). Writing 

1 



F(A6l) 



, „ sin-lT4^(r)i?cos(A6'/2) — 
7^^n^ dr 



(1 - W^{r))-^/^Rcos{A9/2)—W{r) 

dr 



(58a) 
(58b) 

(59) 
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FIG. 3: Correlation Functions. (a,b) Short range correlations, (a) Charge correlation function C(ar) and (b) Zero crossing 
correlation function V{ar) for short range Gaussian field correlator W{ar) — exp{—a?r'^). (c,d) Long range correlations, (c) 
Charge correlation function C(ar) and (d) Zero crossing correlation function T{ar) for long range field correlator W{ar) = 
Jo{ar). Short range correlations: C{ar) in (a) is everywhere negative, indicating that on average positive (negative) charges are 
surrounded by a tight cloud of excess negative (positive) charge, leading to local charge neutralization - this is the characteristic 
field structure of short range screening, see Fig. 1(a). Similarly, T{ar) in (b) is also everywhere negative, indicating that a 
positive (negative) zero crossing on the boundary is likely to be followed by an excess of negative (positive) zero crossings 
that quickly cancel its contribution to the net charge. This structure leads to a linear increase of the charge variance with 
boundary length. Long range correlations: Although for, say, a central positive charge the first negative dip in C{ar) in (c) 
integrates to an excess negative charge of —0.76, indicating substantial, albeit incomplete, local charge neutralization, much 
of this is offset by the following positive peak. As can be seen, the net average field structure consists of rings whose weak 
excess charge alternates positive/negative, with each successive negative ring having marginally greater area, and therefore 
marginally greater negative charge, than the positive ring that follows it. These oscillating cancellations give rise to a slow 
1/r approach to complete screening that is the hallmark of long range screening. This weak structure is too diffuse to be seen 
in Fig. 1(b). For the zero crossing correlation function T{ar) in (d), a positive (negative) zero crossing on the boundary is 
likely to be followed by a long string of zero crossings with a weak excess of positive (negative) signs that reflect the bulk ring 
structure (the boundary runs through one of these rings). This structure leads to an RlnR growth of the charge variance. 



and using the fact that < A6' < tt, so that cos (A6'/2) > = ^47?^ _ r'^/{2R), we have 



rioTrR 1 



2R 



dr 



2 TT 7o Rcos{A9/2) 
d 



Rcos{Ae/2) sin-^ W{r) 



{1 - W^{r))-^/^Rcos{A9/2) 



dW{r) 



dr 

tIoTtR 1 '■^^ 



dr 



2n 



\ dr sin" ^W(r) 



d_ 
dr 



\J\R^ - r2 



W'{r) 
- W^{r) 



(60a) 
(60b) 



Q.E.D. 
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B. Approximate Equivalences 



It was suggested in [9] that the zero crossings of TZ or of I need not necessarily be counted over the actual perimeter 
P of A, but rather that in an isotropic system one could perform the count over any straight line with the same 
length as P. This suggestion, supported by numerical analysis and computer simulations, was based on the notion 
that for short range correlations the geometry of the line - rectangular, circular, straight, etc. - was unimportant, 
as long as over most of its length its radius of curvature was large compared to the screening length A^. Here we 
demonstrate the validity of this suggestion analytically for short range correlations, and examine its validity for long 
range correlations. 

We define an error parameter £ by 



((AiV)')L 
£ 



£1 + 



(61a) 
(61b) 



where 



£1 



no 
' 2 



drrT(r), 



(62) 



p2tvR 

£2 — fioirR I drT{r), 
J2R 



(63) 



£3^ — drsin"! W(r)— ( 
27r Jq dr \ 



2R 



W'{r) 



y/1 - W^{r) 



(64) 



Short range correlations 



For large i?, when W is short ranged r(r) decays to zero well before the upper limit of the integral on the R.H.S. of 
Eq. (62 1 is reached. Extending this limit to 00, it is easily seen that £1 is a constant that does not grow with R, and 
can therefore be neglected; for the same reasons, £2 in Eq. (63) can similarly be neglected. Expanding the square 
root in the R.H.S. of Eq. (64 1, it is easily seen that also £3 does not grow with i?, so it too can be neglected. Thus 

are equivalent in the limit of large R. 
Returning to Eq. (37 1, we have to leading order in r/R 



for short range correlations the linear F formula and the exact result for {Q^ 
This equivalence can be demonstrated more directly. 



R r°^' W'jr) d 
~ TT io ^ y/1 - W^{r) dr 

"""^ 1 + 2 / drT (r) 



arcsin {W{r)) 



(65a) 
(65b) 
(65c) 



(52) 



which is also the leading, i.e. R dependent, term in the linear F formula, Eq 

We therefore conclude, in full accord with [S], that for short range screening t 
used to measure (Q^). 



le linear F formula can, indeed, be 



2. Long range correlations 



For long range correlations we need the asymptotic form of F(r), which for Jo is 

a 1 + sin(2ar) 



F(r) 



(66) 



To leading order the oscillating sin(2ar) /r term makes no contribution, and we have: 
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for £i, 



2ttR 



drr 



27r3r 



for £2, 



So 



aR 



R 
aR 



2ttR 
2R 
InTT, 



dr- 



and for £3 







2R. 



1 
ai? 

— In2-1 . 







-r^ -2R 



2R 



where in evaluating Eq. ( 69a I we use the asymptotic form 

W{r)W"ir) « - 



2a 



Assembling the pieces yields 



iR 



■In 



sin^ (ar) 



M58aR. 



(67a) 
(67b) 

(68a) 
(68b) 

(69a) 
(69b) 
(69c) 

(70) 

(71) 



Thus, unlike the case of short range screening, for long range screening the linear T formula is always larger than 
the true result; the percentage error, however, is not severe, being 8.8% for aR — 10, 6.3% for aR — 100, and 4.9% 
for aR = 1000. 



(^{AN) J exceeds (Q^) because setting L — 2ttR in the integrals over T {Ax) in Eq. (52 1 implies that the 

correlations between zero crossings act along the arc of the circle. But, as already noted in connection with the 
circular T formula, these correlations act only along the chords of the circle. However, except for the diameter, all 
chords are shorter than their arcs, so that L = 2-kR necessarily overestimates the correlation effects. 

W e no w inquire as to what value of L in the upper limits of the correlation integrals appearing on the R.H.S. of 

Eq. (52 1 establishes equality between 1{AN)^\ and (Q^). Writing L — 2pR, we have 



((AiV)^^^ 



= {Q^) + faR, 



f = 



ln(p/2) + 1 - p/tt 
7r2 ■ 



Setting f — and solving for p we obtain 



-ttLW 



2 

7re 



1.01701687... 



(72a) 
(72b) 

(73) 



where LW is the LambertW function [20], and e = 2.718... is the base of the natural logarithm. Eq. (73 1 suggests 
that as a physically attractive approximation L should be set equal to the longest chord, 2R, rather than to the 
longest arc, 27ri?. Indeed, inserting p = I into / yields the excellent approximation £ = — 0.001 16ai? - a forty-fold 
improvement over Eq. (71 1. 
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40 60 

aR 



80 100 



FIG. 4: Practical realization of Jq. (a) Correlation functions. Solid curve W{ar) = Jo{ar), small circles W{ar) equal to the 
practical realization of Jo, Eq. ( |74a[ ); in both cases a = 1 and b — 0.01. The agreement between the two forms is excellent 
out to ar = 100. (b) (Q^) vs. aR. Solid curve the exact result, Eq. (|37fc with Jo correlation function, small circles Eq. (|37fc 
with practical realization of Jo. The inset shows the percentage error obtained using the practical realization. This error is 
less than 0.5% for aR = 50, and less than 1.5% for aR — 100. 



IV. PRACTICAL REALIZATION OF J 







As discussed in Section lie, the long range correlation function Wjg (r) = Jq (ar) arises from a source S (u) that 
takes the form of a ring of zero width, Eq. (36). Here we consider a ring of mean radius a and finite width 6 as a 
practical reahzation, Wpjg (r), of Wj^ (r). 



We obtain Wpj„ (r) from the Van Cittert-Zernike theorem [5T] as 



Wpjo (r) 



1 

abr 
Jo (ar) 



[{a + b/2) Ji {ar + br/2) - (a - 5/2) Ji {ar - br/2)] , 



{ar) 



-(- 

6\2a 



{ar) Ji {ar) + O 




(74a) 
(74b) 



Wpjg (r) closely matches Wj^ (r) over the region < ar < a/b, see Fig. Qa). We can therefore expect for < ar < 
ha/b, and in fact do obtain, see Figs. j4|b,c), a close match between the result obtained using Wpjg in Eqs. (jT]) and 
^ and integrating Eq. ( 37 1 numerically, and that obtained using Wjg . 

A ratio b/a = 100 appears to be experimentally practical, so that the long range screening results presented here 
should be amenable to experimentation. 



V. DISCUSSION 



Smoothed boundaries 



What does smoothing the boundary of the wavefield mean both physically and mathematically? Close examination 
of the derivation of Eq. ( 20a I reveals that all boundary smoothing functions act as filters that smoothly reduce to 



zero not the field amplitude as one recedes from the central point, but rather the number density of singularities. 
How can such a fall-off be achieved? 

A physically meaningful possibility is to devise a source function S that leads to a wavefield with an intrinsic 
singularity density that falls smoothly to zero from some central point. However, not only is there no suitable source 
function known, even in principle, but also such a wavefield, if it could be created, would violate a fundamental 
assumption in the derivation of the charge correlation function C (r) and in the calculation of (Q^) that is used in 
[TD] as well as here, namely that the system is stationary, i.e. that on average it is the same everywhere, that there is 
no central point. Thus, for such a system the results given in , and our extension of them, need not necessarily 
apply. 

So, if an intrinsic average number density of singularities that is independent of position is needed, then how can 
the fall-off required by the smooth boundary in lOj be obtained? The only possibility that comes to mind is the 
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following: Some arbitrary point, the central point, is surrounded by a small circle and a series of narrow concentric 
rings. Suitable fractions of the singularities in these rings are then discarded to create a singularity density profile 
that approximates, say, the Gaussian in |10j , or any other desired smoothing function. 

Of course, an experimentalist would object that not only is such a procedure entirely arbitrary, but also that no 
good could ever come from throwing away data. And, as shown above, she would be right, because the cost of 
boundary smoothing for short range correlations is a loss of all wavefield information. Similarly, for long range 
correlations the magnitude of the calculated fluctuations change with changes in smoothing function, showing that 
the coefficient S in Eq. (47 1 cannot be interpreted as a screening length, as is done in Eq. (12 1. 



B. Deviations From the Linear Law 



Small R for any correlation function 



Surprisingly, perhaps, for sufficiently small R, (Q^) ~ R^ for all field correlation functions for which W (0) =0; 
this includes not only the Gaussian in Eq. ( 17 1 , Jq in Eq. ( 35 ) , the practical approximation to Jq in Eq. ( 74a I , but 



also every other physically realizable short and long range correlation function. W (0) = follows from expansion 
of the Hankel transform in the Van Cittert-Zernike theorem [21J, and the fact that all physically realizable S are 
bounded. The proof of this universal R^ law for i? < 1/a is as follows: Returning to Eq. (37 1 we have for ar <C 1, 



r2R 

Jo 



dr\/AR^ - r^, 



tjttR'^, 
N. 



(75a) 

(75b) 
(75c) 



The physical meaning of this result is that for sufficiently small A = ttR^ there can no be screening because the 
probability of finding the required negative charge within A is vanishingly small. 



2. Large R for long range correlation functions that decay slower than Jo 
(i) If one had a form for W (r) for which W'{r) decays asymptotically as with f3 < ^, then 



22/3+1 (2-/3) 



R 



2-2/3 



(76) 



(li) For /?= i, (Q^) w RlnR, Eq. (|46al 
(in) For /3 > i, (Q^) w R, Eq. 
The proofs of assertions (i) and (m) are as follows: 



For /? > Eq. (13 1 converges, and therefore Eq. il2h holds. 
For /3 < i we rewrite Eq. (37) as 



J^i = — / drv/4i?2 - r2f}2 (r) 
27r Jo 



^1+^2, 
f-A 



^2 = 



27r 



2R 



dr^/AR^ - r'^n'^ (r) , 



(77a) 
(77b) 

(77c) 



where Q. is given in Eq. (8), and 1 ^ aA ^ aR. For Ti we have Ti ~ ^ drO? (r), and since this is subdominant, 
depending only linearly on i?, we discard it without further discussion. Inserting the asymptotic form fi^ (r) w l/r^^ 
into J^2i and writing s — r / (2R), we have 



^2 



i?2-2^ 

22/3-1^ 



(78) 



A/(2fl) 



Passing to the hmit A/ (2i?) ^ we obtain Eq. (76 1 
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3. Physically realizable fields 



How small can /3 be in practice? The answer is that asymptotically, i.e. for sufficiently large R, /3 > 3/2. This 
is demonstrated below using the Van Cittert-Zernike theorem [21] and the fact that all physically realizable sources 
S (u) must be: (i) nonsingular, a requirement that excludes, inter alia, delta functions; (ii) positive definite; and (in) 
strictly bounded, i.e. S {u > Umax) — 0. 

In physically realizable sources whose finite extent is defined by a mask where the intensity falls discontinuously 
from some non-negligible value to zero at the mask edge u — Umax i we have for large r 



W{r) 



uS{u)Jo{ru) du ' 



{Umax) ' S{Umax)cOs{rU 

max 

+ 7r/4)/r / , 



(79) 



so that both W{r) and W'{r) decay as r'^/^, i.e. /3 = 3/2. 

If S{umax) is extremely small the r~^/^ decay of W sets in at values of r that may be so large as to be unmeasurable 
in practice. In that case, in practice one observes /3 > 4. To show this we set Umax equal to infinity and write 



W{r) = W1+W2, 
Wi 



Wo = 



uS{u)Jo(ru) du, 
uS{u)jQ{ru) du, 



where 1 A 1/r. In Wiwe expand S{u) « S{0) + uS'{0) and obtain for large r 

^Ai/2cos(rA + 7r/4) 3 /2'sin(rA + 7r/4) 



Wi 



3(0) 
+5'(0) 



r.3/2 



8 V TT Ai/2r.5/2 



2 A3/2cos(rA + 7r/4) 11 /2 A^/^ sin(rA + 7r/4) 

TT S \ TT ^-5/2 



Using the large argument expansion of Jo in W2 and integrating by parts we have 

sin(ru + 7r/4) 1 cos(ru + 7r/4) 



W2 



2 

71" J A 



uS{u) 



8 (ru)3/2 



du 



y| cos(.A + ./4) [^(o)^^^,(o)]^V2 
sin(rA + tt/4) ( S'(O) + AS"(0) 



p5/2 



2A1/2 



5'(0)A 



1/2 



l sin(.A + ./4) [^(o)+A^^(o)]A-V2 



5(0) 



3(rA + 7r/4)r-3/2Ai/2 + £ sin(rA + 7r/4)r-5/2A-i/2 



-5'(0) 
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cos 



(rA + 7r/4)r-3/2A3/2 _ sin(rA + 7r/4)r-^/2Ai/ 



Summing Wi and yields 



W^(r) « -S"(0)/r3, 



(80a) 
(80b) 

(80c) 



(81) 



(82a) 



(82b) 



(82c) 



(83) 



i.e. (3 = A. If, in addition, S"(0) = 0, the decay of W{r) is faster yet. For example, when all odd derivatives of S{u) 
vanish at the origin, as happens for a Gaussian S {u), W{r) decays exponentially fast. 

Thus, for all physically realizable optical fields (Q^) is bounded at its endpoints by two universal laws: for small 
{Q"^) ^ R^i for large R, (Q^) ^ R- This behavior is illustrated in Fig. 5 for the practical Jq autocorrelation 
function in Section IV. 
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VI. SUMMARY 



Under the assumption of circular Gaussian statistics, the topological charge variance (Q^) of vortices in scalar 
fields [T] and C points in vector fields [5] was analyzed for a circular region of area A and radius R for both short and 
long range correlation functions. 

(z) It was shown, Eq. (10), that when the autocorrelation function W (r) together with the first derivative W' (r) 
go to zero at infinity, Eq. ([l| holds and there is screening. 

(ii) For short range correlations (Q^) grows linearly with R, Eq. (12 1: due to screening, inside A there is local 
charge neutralization, and only partially screened charges nea r th e boundary contribute to the fluctuations. Very 
recent experimental measurements of (Q^) |18j have used Eq. (12 1 to obtain the screening length As, a fundamental 
wavefield parameter. Comparison of with the mean spacing between charges led to the conclusion that the charges 
form small clusters - a conclusion verified by direct imaging of the wavefield. Fig. 1. 

(in) Boundary smoothing [10 was considered for short rang screening, and it was shown that such smoothing does 
not yield useful results, Eq. (32 I. 

(iv) For long range correlation functions (Q^) grows faster than R, and it is not possible to define a screening 
length. For a Jq correlation function (Q^) ^ RlnR, Eq. (46a I. For correlation functions that decay more slowly 
than Jo, - Rp, where 1 < p < 2. 

(v) Although a Jq correlation function is not attainable in practice, it was shown that one can generate an excellent 
approximation valid over an arbitrarily large, but finite range of R, Eq. (74a I; this approximation was shown to yield 
results for (Q^) that are in close agreement with those for Jq, Fig. 4. 

(vi) (Q^) can be calculated using either the charge correlation function C (r), Eqs. ([t]) and (8), or the zero crossing 
correlation function r(r), Eq. (53). An exact calculation showed that, as expected, these two seemingly different 
methods yield the same result for both short and long range screening. 

(vii) For short range screening it was shown that it is also possible to obtain {Q'^) from zero crossing measurements 
made along any straight line of length P ~ 2t:R, rather than only along the circular perimeter of A. For long range 
screening, however, it was shown that this useful simplification no longer holds. 

{via) It was also shown that for every physically realizable wavefield, (Q^) ^ R^ for sufficiently small R, and that 
for sufficiently large R, (Q^) ^ i?. 
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